Dipolar Fermi gases in anisotropie traps 
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The quest for quantum degenerate Fermi gases interacting through the anisotropie and long-range 
dipole-dipole interaction is an exciting and fast developing branch within the cold-atoms research 
program. Recent experimental progress in trapping, cooling, and controlling polar molecules with 
large electric dipole moments has, therefore, motivated much theoretical effort. In a recent letter, 
we have briefly discussed the application of a variational time-dependent Hartree-Fock approach 
to study theoretically both the static and the dynamic properties of such a system in a cylinder- 
symmetric harmonic trap. We focused on the hydrodynamic regime, where collisions assure the 
equilibrium locally. Here, we present a detailed theory, extended to encompass the general case of 
a harmonic trap geometry without any symmetry. After deriving the equations of motion for the 
gas, we explore their static solutions to investigate key properties like the aspect ratios in both real 
and momentum space as well as the stability diagram. We find that, despite the lack of symmetry 
of the trap, the momentum distribution remains cylinder symmetric. The equations of motion are 
then used to study the low-lying hydrodynamic excitations, where, apart from the quadrupole and 
monopole modes, also the radiai quadrupole mode is investigated. Furthermore, we study the time- 
of-fiight dynamics as it represents an important diagnostic tool for quantum gases. We find that 
the real-space aspect ratios are inverted during the expansion, while the one in momentum space 
becomes asymptotically unity. In addition, anisotropie features of the dipole-dipole interaction are 
discussed in detail. These results could be particularly useful for future investigations of strongly 
dipolar heteronuclear polar molecules deep in the quantum degenerate regime. 

PACS numbers: 67.85.-d,67.85.Lm 



I. INTRODUCTION 

Sirice the achievement of Bose-Einstein condensation 
(BEC) with a sample of 52 Cr-atoms jl| the interest in 
dipolar quantum gases has strongly increased @, Q . Re- 
cently, the progress towards quantum degenerate polar 
molecules has pushed this interest even further because 
they possess electric dipole moments of the order of one 
Debye and, therefore, are potential candidates to make 
dipolar effeets more accessible to experiments 0-0] • 

Concerning dipolar bosonic particles, the field is rela- 
tively well understood and has seen a robust development 
with a remarkable quantitative agreement between exper- 
iment and theory. Starting point was the construction of 
a pseudopotential by Yi and You Q. In the Thomas- 
Fermi regime, where the kinetic energy can be neglected 
in comparison with the interaction energy, exact solu- 
tions of the Gross-Pitaevskii equation have been found 
for axially symmetric harmonic traps 0, Ho| . Further 
generalizations to triaxially anisotropie traps even pro- 
vided the first clear experimental signature of the dipole- 
dipole interaction (DDI) in the data for the expansion 



dynamics [llj, |l2|. In the meantime, collisional control 
of chromium has been fully demonstrated by using Fesh- 
bach resonances to increase the relative importance of 
the DDI with respect to the contact interaction. As a re- 
sult, strong dipolar effeets have been observed in BECs 
like the suppression of the characteristic inversion of the 
aspect ratio during the expansion [l3| . In addition, the 
trap configuration could be manipulated to stabilize a 
purely dipolar BEC [l4[ and a new type of 'Bose-nova' 
experiment beautifully revealed a d-wave symmetry in a 
dipolar BEC [Hj]. Besides that, the anisotropie nature of 
the DDI is predicted to shift the Bose-Einstein conden- 
sation temperature in a characteristic way [lfjl fl7j , and 
considcring spinorial degrees of freedom might provide 
an atomic realization of the Einstein-de Haas effect [l8[ . 

In view of fermionic dipolar quantum gases, amazing 
predictions have been made. In the case of homogeneous 
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gases, interesting properties like zero sound 
Berezinskii-Kosterlitz-Thouless phase transition 
nematic phases [IH [23| have been considered, while stud- 
ies of trapped dipolar gases focus on anisotropie super- 
fluidity in three dimensions (2~Ij , fractional quantum Hall 
states [H|, and Wigner crystallization in rotating two- 
dimensional systems [261 ] . 

From the experimental point of view there are differ- 
ent possibilities of realizing dipolar Fermi gases. One 
of them is to use atoms which have large permanent 
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magnetic dipole moments m, such as the isotope 53 Cr 
of chromium, which has a dipole moment of six Bohr 
magnetonsuB and has already been magneto-optically 
trapped [27j or the isotope 173 Yb of ytterbium, which 
has m = 3 [ìb in the 3 P2-state and has already been 
cooled down to quantum degeneracy [28j . In addition, 
recent developments in laser cooling of the isotope 66 Dy 
of dysprosium, with a magnetic dipole moment of the 
order m ~ 10 fiB, promises to increase the variety of 
highly mag netic atomic gases in the quantum degenerate 
regime [29l j . A further exciting possibility is displayed 
by samples of hcteronuclear polar molecules. For them, 
prospects for collisional control through applied electric 
fields 30] indicate that dipolar gases could be explored 
ali the way from the weak- (collisionless) to the strong- 
interaction (hydrodynamic) regime, since this may lead 
to interaction strengths changing by orders of magnitude 
depending on the applied electric field [3l|. This is in 
close analogy to the use of Feshbach resonances to tune 
the contact interaction to unitarity as has been carried 
out with success to observe hydrodynamic behavior in the 
normal phase of atomic Fermi gases [32j . Experimentally 
this is very promising and, recently, 4 x IO 4 fermionic 
40 K 87 Rb molecules with an electric dipole moment of 
about 0.5 Debye have been brought close to quantum de- 
generacy by using stimulated Raman adiabatic passage 
to cfficiently convert the molecules into the rovibrational 
ground state Q. After that, further progress towards 
probing quantum degeneracy has been made, by bring- 
ing this system into the hyperfine ground state as well 
Q. In the meantime, KRb-samples at the temperature 
T = 1.4 Tp, where Tp is the Fermi temperature, have be- 
come available, in which thermodynamic measurements 
led to observation of the anisotropy characteristic of the 
DDI 0]. 

The first theoretical investigations of such a system 
were carried out under the assumption of a Gaussian 
density profile, which is able to capture some of the 
important features as the aspect ratio, but misses the 
correct weak- interaction, collisionless limit [H, [34| . On 
top of that. other approaches have been pursued includ- 
ing [IH IH, [36[ or neglecting (3tJ the possibility of a 
deformation of the momentum distribution. However, 
quantum degenerate heteronuclear polar molecules pos- 
sess strong dipolar interactions which might lead to a 
collisional regime combined with an anisotropie momen- 
tum distribution. For this reason, we have recently devel- 
oped a complete theory for normal dipolar Fermi gases in 
the hydrodynamic regime in the presence of a cylinder- 
symmetric harmonic trap (38|. In the present paper we 
extend our theory to the general case of a triaxial trap, 
which allows to study important aspeets of the physics 
of dipolar Fermi gases such as the radiai quadrupole ex- 
citation as well as to sort out the anisotropie effeets of 
the DDI. 

In the following we treat one-component fermionic 
dipolar quantum gases semi-analytically and tacitly as- 
sume that the gas is in the hydrodynamic regime. 



By adapting a variational time-dependent Hartree-Fock 
method, which was orig inally developed to study nuclear 
hydrodynamics [H, HJ, we are able to obtain a complete 
description of strongly interacting normal dipolar Fermi 
gases which encompasses their static as well as dynamic 
properties. The paper is organized as follows. In Section 
Hi] we outline the variational formalism of hydrodynamics 
applied throughout the paper. In Section Hill we derive 
the action governing the dynamics of the system in the 
case of three different trapping frequencies and extremize 
it with respect to the widths in spatial and momentum 
distributions, obtaining, thus, the corresponding equa- 
tions of motion. Section IIVI provides the dimensionless 
variables which make the physical interpretation of the 
results more enlightening. Then, in Section [Vi we de- 
rive the equilibrium properties such as the momentum 
and real space aspect ratios as functions of the dipolar 
strength and the trap anisotropies. Section fVTl is devoted 
to the low-lying excitations, where we study the oscilla- 
tions around the equilibrium. Following that, we address 
in Section [Vili time-of-flight experiments which represent 
another fundamental issue in cold atoms physics. In Sec- 
tion I Villi we present the conclusion where the main re- 
sults are summarized and further studies of this system 
are discussed. 



II. HYDRODYNAMIC FORMULATION OF 
HARTREE-FOCK THEORY 

Consider a gas containing N harmonically trapped 
fermionic particles of mass M possessing either electric 
or magnetic dipole moments, which are polarized in z- 
direction. The Hamilton operator of such a quantum 
many-particle system is given by 



H 



N r 

E 



2M 



1 



JY 
1 iftj 



(1) 



The first term represents the kinetic energy, which turns 
out to be negligible for Bose systems in the Thomas- 
Fermi regime but is important for Fermi systems, since 
it provides stability against collapse. In general, the trap- 
ping potcntial is harmonic and reads 



tftr(x) 



ÀI 



i 2 2 i 2 2 i 2 2\ 

(uj x x +u y y +oj z z ) 



(2) 



whcr 



and uj z are different trapping frequencies. 



The last term in Eq. ((T|) takes care of the interaction and 
V int (x, — Xj ) denotes the two-body interaction potential 
between particles located at Xj and Xj. 

At very low temperatures, the Pauli exclusion principle 
prohibits s-wave scattering for identical Fermi particles 
and long-range interactions, such as the dipole-dipole in- 
teraction 



^dd(x) = 



a 



dd 



4ir x : 



1 - 3- 



(3) 
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become important. In the case of magnetic dipoles m the 
DDI is characterized by Cdd = Po m2 j with /i being the 
magnetic permeability in vacuum, whereas for electric 
moments we have Cdd = 47rd 2 with the electric dipole 
moment d expressed in Debyes. For most Fermi gases 
which have been realized so far, Vdd is rather irrelevant 
and polarization leads to a degenerate non-interacting 
gas. For polar molecules, however, this is no longer 
valid: moderate electric fields induce dipole moments 
which render the DDI (J3]) a prominent contribution to 
the Hamiltonian in Eq. ([T]). 



A. Center-of-mass expansion 

In the usuai formulation, hydrodynamic studies of de- 
generate Fermi gases are based on closed equations for 
the particle density p(x, t) and the velocity field v(x,t). 
The dynamic properties of this system are determined by 
solving the continuity equation and the Euler equation. 
This set of coupled equations can be obtained from the 
Boltzmann equation for the phase-space distribution (4lj 
or by expanding the equation of motion for the one-body 
density matrix around the center of mass [42j . Let us 
explore further the last possibility in order to illustrate 
some aspects of the method which we apply in the fol- 
lowing. 

Consider the action governing an iV-fermion system 



lo 



A ; l l' ih^ II >l' . 



(4) 



where l^) denotes a Slater determinant built out of one- 
particle orbitals ìpi(x,t) denoted by 



*(x lr -- ,x*r;t) = SD[^(x,t)], 



(5) 



with the energy-level index i taking the values 1 < i < N. 

By extremizing (j4|) with respect to the functions 
ipi(x, t) and !/j*(x, t), one obtains the Hartree-Fock equa- 
tions for the one- particle orbitals ip*(x,t) and éAx,t). 
Combining them yields the equation of motion [43| 



ih 



9/o(x,x'; t) 



2M 



(V x 2 -V x 2 ,)+f/tr(x)-t/ tr (x') 



p(x,x'; t) 



+ [r D (x,i)-r D (x',i)] p(x,x';*) 
Jd 3 r [r E (x,r;i)p(r,x';i)-r E (r,x';i)p(x,r;0] (6) 



for the one-body density matrix 

N 



p(x,x';t) = JJ d Xj**(x'x 2 ,- • •,x A r;t)*(x,x 2 ,- • ■,x. If ;t), 

N 

= £Vi(x,*M(x,t). (7) 



Here, the direct Hartree term, to which only the diagonal 
density matrix contributcs, reads 



r D (x,t) = j d 3 r^ nt (r,x)p(r,0, 



(8) 



while the Fock exchange term, which is given by 

r E (x,x';t) = -^ nt (x,x')p(x,x';t), (9) 

also involves off-diagonal elements of the density matrix. 

In order to obtain the conservation laws corresponding 
to the hydrodynamic equations, we perform an expansion 
around the center-of-mass coordinate X = (x + x')/2 in 
powers of the relative coordinate s = x — x'. In zeroth 
order in s, we obtain from ([6]) the continuity equation 



dt 



•V-j(x,t) =0, 



(10) 



with the particle density p(x, t) = p(x, x; t) and the cur- 
rent density 



j(x,i) 



2Mi 



(V x -V x <)p(x,x';t) 



(11) 



The first order in s yields from ([5]) the Euler equation 
dji(x,t) 



AI 



in 



= -v Xj n°-(x,t)- P (x,t)v Xl c/(x) 



- i o(x,t)V x T D (x,t)+y , d 3 .TV(x,x';t)p(x',x;i)V x K a t(x,x'), 

(12) 

with the non-interacting kinetic stress tensor 



n?.( x ,t) 



fr 2 (Vx-VxOi (Vx-VxQ 3 

' M 2 2 



p(x, x';t) 



(13) 

Introducing the velocity field v(x, t) — j(x, t)/p(x, t) and 
assuming that the trapping potential U (x) is sufficicntly 
smooth, the kinetic stress tensor takes the form [42j 

n°-(x,É) =%P°(x,i)+Mp(x,t)u i (x,t)« i (x,i), (14) 

where the pressure P° (x, t) obeys some equation of state 

P°(x,t) = F(p(x,t)). i 

In case of an irrotational flow, where the circulation of 
the velocity field vanishes due to V x v = 0, the Euler 
equation (|12p can be rewritten in the form 



M 



dv(x,t) 
dt 



p(x,t) 



dp 



3 ,p(x,x';^V(x',x;i) 

L X -. r 



'^ + c/ (x ) + r D (x,t) 



V x M„t(x,x') (15) 



with the transport derivative d/dt = d/dt + v • V. Now 
the effect of the exchange term of the non-local inter- 
action potential ([3]) becomes clear: it breaks the con- 
servation of the circulation of the velocity field v(x, t) 
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and Kelvin's theorem does not hold although we con- 
sider an irrotational fìow. This obvious contradiction is 
a consequence of the fact that it is a priori not possible 
to describe the exchange correlations in terms of dcnsity 
fluctuations alone, Le., fluctuations of the diagonal part 
of the one-particle density matrix. Of course, the true 
exchange correlation is a function of the density alone, 
as a consequence of the Kohn theorem |44j |. and circu- 
lation is conserved. Thus, due the presence of the Fock 
exchange term, the hydrodynamic treatment commonly 
used for dipolar BECs [9] cannot immediately be applied 
to degenerate dipolar Fermi gases. For this reason, we 
propose another approach which preserves the infìuence 
of the non-diagonal part of the one-particle density ma- 
trix [IH , yet assures the conservation of the velocity cir- 
culation. 



B. Common-phase approach 

In this section we discuss the variational time-depend- 
ent approach for a general two-particle interaction poten- 
tial, which will lead to a unified formalism for elucidat- 
ing the hydrodynamic properties of normal dipolar Fermi 
gases. 

In order to study the collective motion of the gas, we 
cmploy a cruciai approximation for the one-particle or- 
bitate ^>i(x, £), namely that they ali have the same phase 



With these definitions the action Q reduces to 



^(x,0 = e lM ^/ ?l |^(x,É)|. 



(16) 



This approximation was introduced before in the context 
of nuclear hydrodynamics [39| and is commonly used in 
hydrodynamic studies (see, for instance, Ref. [401). The 
orbitals | (x, t) | are invariant under time reversion and 
are, therefore, called time-even. 

From Eq. (|16l) and the definition of a time-even Slatcr 
determinant "^(xi,--- , x N' : t) — SD [IV^x, £)|] we ob- 
tain 



*(xi,- • -,x N ;t) 



_ if-[x(x 1 ,t) + ---+ X (x JV ,t)] 



*o(xi,- • -,x N ;t). 
Thus, the one-body density matrix ([7]) reduces to 

p(x,x';£) = e **[x<*' t >-*^.*>U(x )X ';t), (18) 

with po( x , x';t) being a time-even one-body density ma- 
trix. At this point it becomes more evident that the 
present method resembles that of the collective coordi- 
nates applied for fermions, as mentioned in chapter 16 of 
Ref. g|. 

Now the current density, defined in Eq. (jlip. becomes 
j(x, £) = po(x, £)Vx(x, £), allowing for the identification 
of x(x, £) as the potential of the velocity field v(x, £). 



A=-M /d£yd 3 x|x(x,£)po(x,£) 

; 

! 

d£(* |ff|*o). 



Po(x,£) r _ , ..,2 
— o — [Vx(x,£)] 



ta 



(19) 



The first two terms concern the dynamical properties of 
the system and will be shown to give rise to the time 
derivatives in the equations of motion. Notice that inte- 
grating the first term by parts shows that the common 
phase x( x , £) can be seen as the momentum conjugate to 
coordinate po( x , £), which represents the particle density. 
The second term describes the energy associated with the 
movement, Le., the flow energy |45J, given by 



-Eflow(£) 



y /dVolx./)[V\(x. /)]-. (20) 



The last term of Eq. ([ì"9]>. i.e., (\&o|#|*o), consists in 
total of three contributions 

(*0\H\* ) = (*o|^ki„|*o) + (*o|^tr|*o) + (*o|ffmt|*o)- 

(21) 

The first one is the expectation value of the kinetic energy 
operator with respect to |\Po) and gives rise to the Fermi 
pressure: 

-h 2 f 

E kin (t) d 3 x(V x -V x • (V x -VxOpo(x,x';£) 

"(22) 

Notice that the total kinetic energy is given by Er ow + 
(^o|^kin|*o)- For simplicity, the kinetic energy in the 
static case, i.e., (^ol-ffkinl^o)) w ih be referred to as Fermi 
pressure or simply kinetic energy. The second term in 
Eq. (|21j) represents the energy of the external trapping 
potential 



E t r(t) = / d 3 a;po(x,£)C/tr(x). 



(23) 



The interaction energy, given by the third term in 
Eq. (l2~Tj) , contains both the direct and the exchange 
mean- field terms (*o|^int|*o) = E D + E E . The direct 
contribution is given by 

E D (t) - ìy'd 3 a ; dVt/ int (x,x';£) j0o (x,x; £)p (x', x'; £), 



(24) 

while the exchange part reads 

E E (t) = ~^Jd 3 xd 3 x'V ìnt (x, x'; £)p (x, x'; t)p (x', x; £). 

(25) 

Of course, if the interaction energy (^oli? intimo) would 
be a functional of the particle density po( x , £) alone, con- 
servation laws corresponding to the continuity equation 
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and the Euler equation could be immediately derived by 
functionally extremizing the action (|19[) with respect to 
the phase x( x , i) an d the density po( x , i), respectively. 
In the present case, however, one has to extremize with 
respect to the full time-even one-body density matrix 
po(x, x';<). It turns out that the continuity equation re- 
mains unchanged 



dpo(x.,t) 
dt 



= -V-[A)(x,i)v(x,t)] 



The corresponding Euler equation reads, formally, 



M dv^t) 
dt 



3 ,5{y \H\V ) 



d 3 x 



(26) 



(27) 



so that the proposed approach is circulation conserving. 



C. Wigner phase space 

In the preceeding section we have derived a set of equa- 
tions which could be applied to study the hydrodynamic 
excitations of a dipolar Fermi gas. Nevertheless, the 
equations ([26]) and (|27|) are not yet closed due to the 
lack of knowledge of the non-diagonal terms of the den- 
sity matrix. Therefore, we resort to a simpler procedure: 
we extremize action (|19[) with respect to appropriate vari- 
ational quantities, which bear information on both the 
diagonal as well as the non-diagonal part of the interac- 
tion. Of course, we do loose information in this process 
because our variational approach may not be as precise 
as the solution of the complicated equations ©, but, on 
the other hand, it gives access to both the sfaticai and 
dynamical properties of dipolar Fermi gases beyond the 
perturbative regime in a quite simple and clear way 

In order to calculate each of the terms in the action 
(Tì~9l , we change to the Wigner representation of the time- 
even one-body density matrix, which is defined according 
to 

i*(x,k;t) = yd 3 S po(x+|,x-|;t) e- 4ks .(28) 

The inverse transformation reads 

p (x,x';t) = J ^^ (^±^ jk ;t)e ik (— ').(29) 

In the Wigner representation ali quantities of interest can 
be expressed in terms of the Wigner function (|28p . For 
example, the particle density is given by 



po(x,t) = p (x,x;t) 



d 3 k 
(2^)3 



H)(x,k;t), (30) 



and the momentum distribution is obtained via 

q3„ 



f d 3 x 

J (&f i/0(X ' k;t) ' (31) 



With the help of these quantities, the kinetic energy ([22]) 
and the trapping (|23j) energy can be written as 



_ , . f d 3 xd 3 k , 
Euin(t) = J 3 ^ o( - x ' k;t ) 



h 2 k 2 



E tl {t) 



(27T) 

d 3 xd 3 fc 



(32) 



(27T) 



2M ' 

., !/ (x,k;t)C/ tr (x), (33) 



respectively. Accordingly, the direct tcrm in (|24|) . which 
represents the mean-field dipolar potential energy, reads 



d 3 xd 3 kd 3 x'd 3 k' 
2(2tt) 6 



v (x, k;i) %a(x- x')i/ (x', k';t) . 

(34) 



This tcrm is dctcrmincd by the particle density (|3"0")) 
alone and was first considered to analyze the equilib- 
rium [33| and the dynamical [3~i| properties of a cylinder- 
symmetric system by adopting an isotropie momentum 
distribution and a Gaussian trial particle density. 

In contrast, the exchange interaction term (|25[) . given 
by 

e; . . fd 3 Xd 3 kd 3 sd 3 k' , v , . , . , v , , . 
EÌM = 7^-6 ^(X,k;t)y dd (s)z/ (X,k';t) 



x e 



2(2tt) 6 

is-(k-k') 



(35) 



is rather linked to the momentum distribution (|31[) and 
vanishes if it is isotropie. The importance of this term 
in dipolar Fermi gases was only recently recognized in 
Ref. [3a], where it was shown to lead to a deformed mo- 
mentum distribution. Following this important investi- 
gation, further effeets have been studied by taking this 
term into account, like its influence on the dynamical 
properties of cylindrically trapped systems in the colli- 
sionless re gim e [3 61 . For homogeneous gases, zero-sound 
dynamics [19i . |2C| as well as quantum phase transitions 
in two dimensions (2l1 | have been investigated by consid- 
ering also this exchange contribution ([35]) . 

In the following we extend the static, semi-classical 
theory of dipolar Fermi gases to a dynamical theory in 
the hydrodynamic regime by including the Fock exchange 
term in a naturai way. Actually, as is clear from its 
derivation, this theory can also be successfully applied 
to fermionic systems with other types of long-range in- 
teractions. 



III. EQUATIONS OF MOTION 

In order to study the trapped dipolar Fermi gas, we 
adopt for the common phase the harmonic ansatz 

X(x,t) = - [a x (t)x 2 +a y (t)y 2 +a z (t)z 2 ] , (36) 

which is able to capture different excitation modes by 
specifying the form of the potential for the particle flow. 
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Furthermore, we use an ansatz for the Wigner function 
which rcscmbles that of a non-interacting Fermi gas in 
the semiclassical approximation. With this we cope with 
the mairi cffcct of the DDI that the gas is stretched in 
the direction of the polarization. This ansatz is a gen- 
eralization of the one presented in Ref. (35j . which has 
the form of the low-temperature limit of the Fermi-Dirac 
distribution 



i*(x,k;f) =61 



E 



E 



K t (ty 



(37) 



with 6 (x) denoting the Heaviside step function. Accord- 
ing to Eq. (|37[) . the parameters Ri and Ki represent the 
largest extension in the i-th direction of the density and 
momentum distribution, respectively. They will, thcrc- 
fore, be called the Thomas-Fermi (TF) radius and Fermi 
momentum in the i-th direction, respectively. 

Now we are in a position to evaluate the action (fH))) 
as a function of the variational parameters. Introducing 
additionally the chemical potential /i as the Lagrangc 
parameter, which is responsible for particle number con- 
servation, the action reads 



f , R K ! M v-^ r . „„ 

ti 1 



\- h 2 K? _ 3 



2M 

t 2 



' 1 TV RJ 1 \K X ' K v 



- / df n(t) I ~~~~ — - V 



48 



(38) 



where • denotes geometrical average and the Constant co 
is given by 



co 



2 10 C, 



3 4 • 5 • 7 • ir 3 



0.0116 a 



dd • 



(39) 



The DDI is refiected in the anisotropy function /(x, y), 
which is defined as (see Fig. [T|) 

/(x,y) = i + ^y g) -f y j ) , (40) 
(1 — y )vl — x z 

where F(ip, q) and E(cp, q) are the elliptic integrals 
of the first and second kind, respectively, with <p = 



arcsin yl 



and q 2 = (1 - y 2 )/(l - x 2 ). This func- 



tion has often appeared in the literature of dipolar Bosc- 
Einstein condensates [TU, [T3|. Notice that f(x,y) is 
bounded between 1 and —2 passing through at x = y = 
1. This refiects the fact that the DDI is both partially 
attractive and partially repulsive, depending on whct lici- 
tile dipoles are head-to-tail or side-by-side to one an- 
other, respectively. For polarization along the symmetry 
axis, therefore, the interaction is dominantly attractive in 
cigar- (prolate) and repulsive in pancake-shaped (oblate) 



f(x,y) 



fs{x) 




X 



FIG. 1: (Color Online) Top: Anisotropy function f(x,y) from 
(|40[) . which is bounded between 1 for small values of either x 
or y and —2 for large values of both arguments. Notice the 
symmetry f(x,y) = f(y,x) and that f(x,y) reduces to fAx) 
(black curve) in the case of cylindrical symmetry [Ì3 , [TtJ. 
Bottom: f s (x) as a function of x. Notice that f 3 (x) changes 
its sign at x = 1. 



systems. We remark that this is not the only possible 
way to define the anisotropie dipolar function [35]. The 
motivation for our choice is twofold. From the mathe- 
matical point of view, much is known about the function 
(j40| due to the extensive research on dipolar BEC's. For 
example, a good amount of tcchnical information can be 
found in the appendix of Ref. [12]. In addition, it is 
physically appealing to state that the anisotropy of har- 
monically trapped, dipolar particles is determined by one 
and the same function no matter if they are of bosonic 
or fermionic nature. 

In the following, we will use the indexes 1 and 2 to 
denote a derivative with respect to the first and second 
argument. Note that, whereas f(x,y) is symmetric with 
respect to exchanging the first and second variables, this 
is not the case for the functions fi(x, y) and fi(x, y). Fur- 
thermore, in the case of x = y, the anisotropy function 
f(x,y) reduces to @,[Ìl,liÌ] 



/.(*) 



1 + 2x 2 - 3x 2 E(x) 



1 



together with the abbreviation 



E(x) 



; tanti" 1 vT 



< x < 1 



(41) 



(42) 



Vx 1 



Inspecting the action 



tan 1 v ' x 2 — 1; x > 1 

one perceives that the 
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anisotropy function occurs twice. The first time as a 
function of the TF radii Ri, due to the direct term, and 
the second time as a function of the Fermi momenta Ki, 
due to the exchange term. Since /(l, 1) vanishes, this 
term only contributes in the case of deformed Fermi sur- 
faces. In other words, the absence of this term would 
lead immediately to a spherical momentum distribution. 

The equations of motion follow from extremizing the 
action (|38[) with respect to ali variational parameters 
ckì,Rì,Kì as well as the Lagrange parameter fj,. The lat- 
ter assures particle number conservatimi 



where the numerical Constant Cd reads 



A 3 A 3 = 48 A 



(43) 



and is explicitly given by 
1 ^ h 2 Kf 21c N 
3< _ 



2M 



R 



K x 



A„ 



After some simple though tedious algebra, one obtains 
the following equations for the Fermi momenta 

h 2 Kl 



2M 

2M 

h 2 K 2 
2M 



-y 

3 ^ 

l 

-E 

3 ^ 

l 



h 2 K 2 


48Ac K z 


2M 


2R x R y R z K x 




48Ac K z 


2M 


2R x R y R z K y 


n 2 Kf 


48Ac K z 



fi 



K-r 



K, 



2M 

l 

48Ac 



A, K_z_ 

K X Ky 



2R x R y R z K x \K 



A, 



A, 



2R x R y R z Ky 12 \K X ' K y 



(45) 



These equations are clearly redundant, so we drop the 
third of them. Together with the condition for particle 
conservation (|43|) . they determine the Fermi momenta 
Ki as functions of the TF radii Ri, so that we have three 
independent equations to solve for three variables. 

The equations of motion for the variational parameters 
ai are simply given by 



2f 



Cd 



• 5 • 7 • 7T 2 



0.2791. 



(49) 



The first term on the right-hand side of equations (|4"?| 
accounts for the harmonic trap, the second is due to the 
Fermi pressure, and the third represents the DDI contri- 
bution, which will be discussed in more detail in the next 
section. 

Having collected the equations of motion for ali the 
variables, we can attempt to interpret Eqs. (|4"3"|) . (|4"5|) . 
and (|47|) physically. In the case of a spherically sym- 
metric momentum distribution one could neglect the ex- 
change term and set to zero ali terms which involve 
f(K z /K x ,K z /K y ) and its derivatives. Thus, we could 
solve (|4"3")) and (ff51 for the Fermi momenta and obtain 



V48AT 



K x = K y = K z = Kp = — = 

R 



(50) 



Inserting this result in Eqs. (|4"T|) . we, then, would have a 
set of equations of motion for the TF radii which stem 
from a potential V(R X , R y , R z ). The problem would be 
reduced to study the movement of a fictitious particle 
under the influence of this potential. Due to the pres- 
ence of the Fock term, however, it is not possible to solve 
Eqs. (j43"l) and (|45|) directly, so they have to be solved si- 
multaneously with Eqs. (|4"T|) . Thus, we conclude that the 
exchange term modifies the constraints in an anisotropie 
manner such that one has to give up the notion of an un- 
derlying potential V(R X , R y , R z ). Furthermore, it is the 
presence of the Fock exchange term, as first pointed in 
[35j |. which deforms the Fermi sphere into an ellipsoid in 
the case of a cylinder-symmetric trap. As we shall see 
below, this deformation turns out to remain ellipsoidal 
for triaxial traps and plays an important rolc for detcr- 
mining both equilibrium and dynamical properties of the 
system. 



Ri 
Ri 



(46) 



and are used to derive the equations of motion for the 
TF radii: 



Ri = —u 2 Ri + 



^ 3M 2 Ri 

3 



48Ac 



Qt (R,K) 



(47) 



Here the auxiliary functions are given by 

Cd 



Qx (r,k) 

Q v (r,k) 
Qz (r,k) 



x 2 yz 

Cd 



xy 2 z 

Cd 



xyz< 



A 



x y 

z z 

x y 

z z 

x y 

z ' z 



/,( 
/,( 



x y 

z z 



x y 

z z 
x y 
z ' z 



k 

K 

\z zJ 



j, | k z k z 

k x k x 



(48) 



IV. DIMENSIONLESS VARIABLES: 
CYLINDRICAL SYMMETRY OF THE 
MOMENTUM DISTRIBUTION 

Bcfore we explore the physical consequences of the 
equations of motion for a trapped dipolar Fermi gas, let 
us briefly discuss the non-interacting case, which will pro- 
vide us with adequate units for the quantities of interest 
throughout this work. 

Denoting the Fermi energy of a non-interacting 
trapped Fermi gas by Ep, its chemical potential takes 
the form 



fi {a) = E F = fEJ{6N)^ 



(51) 



and the Fermi radii and momentum read, respectively, 



i?i 0) = 



' 2E F 
Muj 2 



K F 



2ME F 



h 2 



(52) 
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This motivates to express the TF radii Ri in units of 
R^ and the Fermi momenta Ki in units of Kp- Defining 
Ki = Ki/Kp and Ri = Rì/r[°\ the condition for the 
particle number conservation (I43|) reduces to 



—3^3 

R K 



1. 



(53) 



Thus, the equations of motion for the TF radii in the 
dimensionless notation will be written in terms of the 
ratios 



Ri 
~R~ 



R,. u>i 



(54) 



outlining the role played by the trap frequency ratios 
K = w z /u x and X y = L> z /u> v . 

Before solving the equations of motion, we can aheady 
obtain important information by considering the symmc- 
tries of the total energy in the static case, Le., with the 
velocity potential x se t to zero. As a function of the 
variational parameters, the energy then reads 



E 



NE F 



52 \ 2e d dCd 
Hi I = g — 

R 



f 



Rx Ry ^y 



' \K X , Ky / 



,(55) 



Rz Rz 

where the dimensionless dipolar strength edd is given by 



£dd 



47T v h 5 



(56) 



Notice that e d d depends on the particle number and on 
the trap frequencies. This is in contrast with the dipo- 
lar Bose gas, where, for a system with s-wave scatter- 
ing length a a , the corresponding dimensionless dipolar 
strength e d d is given by 



fdd 



C dd M 
\2-KÌì 2 a s 



(57) 



The difference between (|56|) and (1571) has important con- 
sequences for the behavior and the tunability of the sys- 
tem, as will be explained in more details below. 

As a consequence of the symmetry of the anisotropy 
function f(x,y) = f(y,x), the energy (|55)l possesses the 
same symmetry with respect to the piane xOy in both 
K- and i?-space. On the one hand, this implies that, in a 
cylinder-symmetric trap, where we have X x = X y , the ex- 
trema of the energy satisfy R x = R y . On the other hand, 
since the trap geometry does not innuence the exchange 
contribution to the total energy, we conclude that the 
momentum distribution of a dipolar Fermi gas remains 
cylinder-symmetric even in the case of a triaxial trap ge- 
ometry, Le., one has K x = K y , provided that the dipoles 
are aligned along the Oz direction. Therefore, in the 
expression above, f(K z /K x ,K z /K y ) can be simplified 



to f s (K z /K x ) without loss of generality. Furthermore, 
noticing the limits 

hm x—f(x, y) = lim y—f(x, y) = — H — ¥ — , 

y^x tJx y^x oy 2(1 — x z ) 

(58) 

we conclude that Eqs. (|45p reduce to the single condition 
Kl-kl = e dd c(K,K x ,K^ (59) 



with the function 

c(n,K x ,k z ) = ^§ 

R 

(60) 

which was aheady found in our previous work concerning 
the cylinder-symmetric trap [38j . In that work, we traced 
Eq. (|5"9")l back to the Fock exchange term. 

We emphasize that the cylindric symmetry in momen- 
tum space also holds in the dynamic case, since neither 
modulating the trap frequencies nor turning them off af- 
fects the symmetries of the exchange term. 

The equation (|47|) for the Thomas-Fermi radius in the 
i-th direction can, thus, be written as 



1 d 2 Ri 



k 2 

= -R i + Y / -^--e dd Q i (K,K x ,K z ), (61) 
. óRì 



with the corresponding simplifications in the Qi-func- 
tions ([15]). 

The equations (|53[) . (|59)l . and (|6T|) describe both the 
static and dynamic properties of a triaxially trapped 
dipolar Fermi gas in the hydrodynamic regime and rep- 
resent the main result of the present article. In what fol- 
lows we shall explore their solutions in different cases of 
interest such as the conditions for stable equilibrium, the 
low-lying oscillations around the equilibrium positions, 
and the expansion of the gas after release from the trap, 
Le., the time-of- flight dynamics. 

In terms of the dimensionless quantities introduccd 
above, the chemical potential (|4"4")) is given by 



E F 



E 



K 2 7c d e dd 
3 8R 



fi 



Rx^-x Ry^y 



Rz 



Rz 



fs 



(«2) 

This expression can also be obtained through the virial 
theorem for dipolar gases [33| 

£ ki „ - E tI + 3E int /2 = 0, (63) 

together with the scaling relation for the particle number 

Nfi = 5£ kin /3 + E* + 2E int . (64) 

Now that we have explained in detail how our equa- 
tions of motion arise and how they are expressed in di- 
mensionless units, we are allowed to drop the ~ 's and 
obtain a cleaner notation without any danger of misun- 
derstandings. 
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FIG. 2: (Color Online) Aspect ratio in real space R X X X /R Z 
for a cylinder-symmetric trap with X x = X y = 7, 6, 5, 4, 3, 2, 1 
(top to bottoni). The upper brandi (continuous) corresponds 
to a locai minimum of the total energy, while the lower branch 
(dotted) represents an extremum but not a minimum. 




FIG. 3: (Color Online) Aspect ratio in real space R X X X /R Z 
for a triaxial trap with X x = 5 for various values of X y = 
3, 4, 5, 6, 7. Notice that, for a fixed \ x , making X y larger cor- 
responds to flattening the trap in the direction perpendicular 
to the dipoles, allowing for stable configurations for larger 
interaction strengths. 



V. STATIC PROPERTIES 

The static properties of a dipolar Fermi gas are ob- 
tained from Eqs. (|61[) by requiring the left-hand side to 
vanish. However this only gives us the conditions for 
an extremal mean-fìeld energy. Since the dipolar inter- 
action also contains an attractive part, it is useful to 
have a criterion for deciding whether a given state, Le., 
a point (R, K x , K z ) in the five-dimensional space of vari- 
ational parameters, is stable or unstable. To that end 
we turn to the total energy, given by Eq. (|55p. which 

3 3 

shall be minimized under the constraint R K — 1 due 
to particle number conservation. A dimensionai analysis 
of the energy (|55p shows that the system cannot have 
a global minimum for any non-vanishing edd- This can 
be roughly seen by noticing that the stabilization comes 
from the factor K 2 ~ Rr 2 whereas the dipolar interac- 
tion goes with iì~ 3 , rendering the energy not bounded 
from below. Nonetheless, for weak enough interactions a 
locai minimum might exist, to which the system would 
return after a small perturbation. The regions satisfying 
this property will be called stable, while infìection points 
and locai maxima will be denoted unstable equilibrimi! 
points. The mathcmatical criterion behind this classifi- 
cation is given by the eigenvalues of the Hessian matrix 
associated with the four effectively independent variables 
of the problem. 

One of the consequences of the unboundedness of the 
internai energy is that, for each value of the interaction 
strength edd for which the system presents a stable con- 
figuration, there is also another unstable one. This can be 
seen by considering the aspect ratio of the cloud, which 
is depicted in Fig.[2]for diffcrcnt values of the trap aspect 
ratio X x = X y = A as functions of edd- Here, we recog- 
nize that the stable branch (continuous) of the real space 
aspect ratio starts at edd = with R x = R z = 1 and 
extends itself until the value ejj"*, where it meets the un- 
stable branch (dotted). For edd > Cdd*' 110 stationary so- 



lution for the equations (|61[) exists. The unstable branch, 
on the other hand, possesses a vanishing aspect ratio for 
£dd = 0. This is due to the fact that the DDI tends to 
stretch the sample along the polarization direction. For 
a small value of edd, the unbounded energy solution is 
obtained with R x — > and, consequently, R x /R z — > 0, 
although the TF radius in the axial direction R z remains 
finite. We remark that the upper branch corresponds 
to a locai minimum of the energy such that the Hessian 
matrix has only positive eigenvalues, while the lower one 
is an extremum but not a minimum, corresponding to a 
Hessian matrix with at least one negative eigenvalue. The 
corresponding graph for a dipolar BECs in the Thomas- 
Fermi regime bears a cruciai difference: unstable Solu- 
tions only become available for edd > 1 @ ■ The physical 
reason for this effect is that in dipolar BECs the stabi- 
lization comes from the contact interaction, which scales 
with R~ 3 , just like the DDI. 

In order to study the effect of a triaxial trap on the 
static properties of a dipolar Fermi gas, we explore fur- 
ther the symmetry f(x,y) — f{y,x) of the anisotropy 
function as defined by Eq. (|4T)]) . Due to this symmetry, 
we only need to discuss the aspect ratio R X X X /R Z since 
the properties of R y X y /R z can be obtained by analogy. 
As indicated in Fig. |3l varying X y for fixed A^ clearly 
affects the stability of the system. For X y > X x , stable 
solutions are admitted for larger values of edd, i-e., e^ 1 ' 
is shifted to the right, whereas in the case X y < X x , 
decreases. This reflects the fact that more oblate traps 
tend to allow for larger edd because they favor the re- 
pelling part of the interaction. Another worth remarking 
feature in Fig.[3]is that reducing X y for fixed A^ reduces 
the value of e d d' much more than it is enlarged by in- 
creasing X y . 

Concerning the aspect ratio in momentum space, we 
have studied its dependence on the dipolar strength edd 
and found an analogous behavior to the one in real space. 
This goes back to the property of the function f s (x) of 
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FIG. 4: (Color Online) Aspect ratio in momentum space 
K x /K z for X x = 5 and X y = 3,4,5,6,7. The upper branch 
(continuous) corresponds to a locai minimum of the total en- 
ergy, while the lower branch (dashed) represents an extremum 
but not a minimum. 




Stable 



0.1 0.5 1.0 5.0 10.0 50.0 100.0 

FIG. 5: (Color Online) Stability diagram of a dipolar 
Fermi gas. The middle, black curve represents the cylinder- 
symmetric case X y = X x , while the upper, red one represents 
the case X y = 5X X and the lower, blue curve is for X y — X x /5. 



othcr. 



changing sign at x = 1, so that the minus sign in front 
of f s (K z /K x ) in the total energy partially compensates 
its dependence on the reciprocai momentum aspect ratio 
K z /K Xl and its behavior with respect to edd turns out 
to be analogous to one in real space. This is explicitly 
shown in Fig. [4j where the aspect ratio in momentum 
space K x /K z is plotted as a function of edd for X x = 5 
and X v = 3,4,5,6, and 7. The main difference, which 
appears in momentum space, is the observation that the 
unstable solution converges to a finite value of the as- 
pect ratio as the interaction strength edd approaches zero. 
This reflects the fact that the collapse is a real-space 
phenomenum which is dominated by the shrinking of the 
radiai Thomas-Fermi radius R x , while the axial Thomas- 
Fermi radius R z remains finite. As the momentum-space 
variables are accounted for only by the constraint (|53|) 
and the condition for momentum deformation (|59|) , both 
K x and K z diverge as edd approaches zero in the unstable 
branch, but their ratio always is finite. 

To conclude our investigation of the static proper- 
ties of a trapped dipolar Fermi gas, we have calculated 
the stability diagram for the cylinder-symmetric case 
A x = Aj,, where we obtain similar quantitative results 
as in Ref. (35|, and for X x ^ X y , where the lack of axial 
symmetry has a considerable influence. The results are 
presented in a log-log plot in Fig. [5] If we consider a 
situatimi in which X y = 5X X (red, upper curve), we do 
not obtain a large variation with respect to the cylinder- 
symmetric case Aj, = A^ (black, middle curve). On the 
contrary, if we take Aj, = A^/5, appreciable differences 
can be noticed as X x increases. This can be understood 
if one realizes that it is the weaker trap frequency which 
determines the highest value of edd , for which the system 
remains stable. Therefore, by enlarging X y with respect 
to X x one obtains a smaller difference with respect to the 
case X y = X x than by reducing it, explaining the effect 
already anticipated in Fig. [3] Also remar kable is the fact 
that, for small A's, the three curves lie very close to each 



VI. 



LOW-LYING EXCITATIONS 



The low-lying excitations of a dipolar Fermi gas are 
studied in this section by linearizing the equations of mo- 
tion around the equilibrium. This is done by assuming 
that at time t the following ansatz is valid 



Rt(t)=Ri(p)+rne 



;nt. 



K t (t)=K l (Q) + ( z e lUt , (65) 



where r\i and Q denote the small amplitudcs in real and 
momentum space, respectively, while fi represents the 
frequency of the oscillations. Due to the cylinder symme- 
try in momentum space, there are only two independent 
momentum-space amplitudes Ce and £ 2 , while three inde- 
pendent real-space amplitudes occur for a general three- 
dimcnsional trap. 

Bcfore we can derive a matrix equation for the ampli- 
tudes in real space, we must obtain the C's as functions 
of the rji's. To that end, we expand Eqs. (|5"3"|) and (|5T))) 
up to first order in the amplitudes and get 



(66) 



with the abbreviations 

K 1 + K 2 Z - e dd K z C. Kz 



W x 

w z 



AK* + 2K X - 3e dd K z C K , 

2K 2 Z - e dd K z C, Kz 
'4K* + 2KÌ-3e dd K z C Kz 



(67) 



To make the notation more succint, we have introduccd 
here the shorthand A t K z = 9A(R,K X ,K Z ) /dK z to de- 
note a partial derivative of the quantity A (R, K x , K z ) 
with respect to K z evaluated at equilibrium. These re- 
sults show that the presence of the dipolar exchange 
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Kx 
K. 



FIG. 6: (Color Online) The lower (red) curve shows the ratio 
of the amplitudes C^/Cz as a function of edd for \ x = \ v = 7. 
For comparison, the stable branch of the equilibrium aspect 
ratio in momentum space against edd for X x = X y — 7 is 
depicted by the upper (blue) curve. 



term drives the momentum oscillations anisotropie. In 
Fig. [6] we plot the ratio Gx/Gz a s a function of edd for 
X x = X y = 7. We show also the stable branch of the cor- 
responding equilibrium aspect ratio in momentum space, 
represented by the upper (blue) curve. In order to appre- 
ciate the meaning of this curve, let us consider a typical 
experimental situation with N ss 4 • IO 4 KRb molecules 
and trap frequencies of (uj x ,u) y ,uo z ) = 2ir (40,40,280) 
Hz. By using an external electric field and tuning the 
electric dipole moment to d = 0.2 Debye, one obtains 
the dipole-interaction strength edd ~ 0.43 which leads to 
an oscillation anisotropy of Gx/ Gz ~ 0.87 and an equilib- 
rium momentum deformation of K x /K z ss 0.93. More 
striking effeets result for a stronger interaction. Con- 
sidering an electric dipole moment of d = 0.57 Debye 
yields a larger value edd ~ 3.53 and, therefore, also 
larger anisotropies for both the momentum oscillation 
Gx/Gz ~ 0.45 and equilibrium momentum distribution 
K x /K z ss 0.64. These results exhibit clcarly the effeets 
of the exchange term on the low-lying oscillations and 
make room for a clear detection of the DDI in ultracold 
degenerate Fermi gases. 

Linearizing the equations for the TF radii (|61[) we ob- 
tain 



if_ 2K 2 

2 ~~ 



K 2 



3Rf 



E 



p 



RiRj 



where we have introdueed the shorthand 



(68) 



P 



2Kt_ 
3 



Kt + AK 2 K 2 - (2K 2 + K 2 ) e AA K z C, Kz 
AKl + 2K1 - 3e dd K z C, Kz 



(69) 

which approaches the value 2/3, as edd goes to zero. 

With this the study of the low-lying oscillations in a 
dipolar Fermi gas has been reduced to the eigenvalue 
problem (|68l) : the oscillation frequencies i7 are given by 
the square root of the corresponding eigenvalues and the 
eigenmodes describe the real-space motion during the os- 
cillations. 



For a non-interacting Fermi gas, this formalism recov- 
ers the oscillation frequencies of a triaxial trap as the 
solutions of the equation 

3fi(°) 6 - 8^°> 4 (^ + <4 + + i^ (0 %W v + <4<£+ 

uj 2 y z ) - 48uluy, - 0. (70) 

This result is in agreement with Ref. 14 7|. where a 
deeper analysis, initially devised for BECs [4g], is car- 
ried out. It is shown there that, despite the lack of an 
obvious spatial symmetry, the wave equation for the hy- 
drodynamic modes is separable in elliptical coordinates. 
We remark that the solutions of Eq. ([7D1 reduce to the 
respective frequencies in the presence of cylindrical [42| 
or spherical symmetries [41| , where this problem was first 
tackled. 

In the following, we discuss separately the effeets of 
the DDI in cylindric and triaxial traps. The modifica- 
tion introdueed in Eq. (1701 due to the inclusion of the 
DDI makes this equation too cumbersome to be displayed 
here. The same is true for the corresponding solutions. 
For this reason, we shall provide detailed expressions for 
the three oscillation frequencies only in the case of cylin- 
der symmetry. 



A. Oscillation frequencies in cylinder-symmetric 
traps 

In the presence of cylinder symmetry, we find three well 
characterized oscillation modes: one two-dimensional 
mode, the radiai quadrupole, and two three-dimensional 
ones, the monopole and the quadrupole modes. 

The first mode we consider is the radiai quadrupole 
mode, depicted in Fig. [7^,). It is characterized by a van- 
ishing amplitude in the Oz-direction, while the oscilla- 
tions in Ox- and Oy-directions have the same amplitude 
but are completely out-of-phase. We find the frequency 
fì rq to be given by 



firq = uj x <2 + 



3cd£dd R X X 2 

RÌRr RI 



2 (R 2 Z - R 2 X X 2 ) - {AR 2 Z + R 2 X \ 2 ) f s (R X X/R Z 
(R 2 - R 2 X \ 2 ) 2 



1/2 



(71) 



where the TF radii R x and R z correspond to the static 
values calculated in Section IVI 

The radiai quadrupole mode can be experimcntally ex- 
cited by adiabatically deforming the circular trap in the 
xOy-pl&rie into an ellipse and suddenly switching off the 
deformation. In the case of a Fermi gas with contact in- 
teraction, the radiai quadrupole mode was used to probe 
the transitimi from the collisionless to the hydrodynamic 
regime throughout the BEC-BECS-crossover [49j . For a 
dipolar Fermi g RS, ci similar experiment could be thought 
of, where controlling the collisions through applied elec- 
tric fields would play the role of a Feshbach resonance to 
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FIG. 7: (Color Online) Oscillation modes in the cylinder- 
symmetric configuration for A = 5. The oscillation frequencies 
of the radiai quadrupole Q rq (black), the monopole Q+ (blue, 
dark grey), and the three-dimensional quadrupole mode f2_ 
(red, light grey) are plotted as functions of the dipolar in- 
teraction strength edd in units of u x . The inset shows the 
behavior of the corresponding eigenmodes. Part a) character- 
izes the radiai quadrupole eigenmode, while b) and c) refer to 
the monopole and quadrupole mode, respectively. 
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FIG. 8: (Color Online) Low-lying oscillation frequencies in 
units of the corresponding non-interacting values as func- 
tions of A for edd = 0.8 (dotted) and edd = 1.2 (continuous). 
The monopole and quadrupole modes are displayed in blue 
(dark grey) and red (light grey), respectively, while the radiai 
quadrupole mode is shown in black. 
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tune the system ali the way from a ballistic to a hydro- 
dynamic behavior in the normal phase. Fig. [7J depicts 
O rq as a function of edd for A = 5 in units of u) x . We find 
that, for a given A, the frequency fi rq is quite insensitive 
to changes in the interaction over the range of values in 
which the gas is stable. Although we have varied the trap 
anisotropy A from 0.2 up to 20, no significant alteration 
of this behavior could be detected. Fig. [5] shows the de- 
pendence of f2 rq on A for edd = 0.8 and edd = 1.2 in units 

of its non-interacting value, Le., Q$ = y/2~L) x , which is 
directly given in Eq. (|TTj) by setting edd = 0. We re- 
mark that the function in the second line of Eq. (|71l) is a 
function of the ratio R x \/ R z alone which approaches the 
value —16/35 as R X X/R Z tends to 1, so that no divergence 
arises for R z = R X X. We would also like to point out 
that, despite the fact that the radiai quadrupole mode is 
inherent to cylinder-symmetric systems, its calculation 
rcquires that one starts from a triaxial framework, which 
is then specialized to axial symmetry. This is the rea- 
son, why this important mode was not explored in initial 
studies of dipolar Fer mi g ases in the hydrodynamic [38[ 
or collisionless regime [36j. 

We now concentrate on the three-dimensional mono- 
pole and quadrupole modes. The first, also known as 
breathing mode, is a compression mode characterized by 
an in-phase oscillation in ali three directions and is de- 
noted with an index +. The second, in analogy with the 
radiai quadrupole mode, is an out-of phase oscillation in 
radiai and axial directions and is denoted with an index 
— . In a spherical trap, these modes are decoupled from 
each other, but, here, they are coupled due to the cylin- 
der symmetry of the trap. Their frequencies are given 
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The dependence of the mono- and quadrupole oscilla- 
tion frequencies on the DDI strength edd for a fixed trap 
anisotropy A is shown in Fig. [7] We find that the frequen- 
cies behave for different values of A > 1 qualitatively like 
in Fig.[jj where we have A = 5: The monopole frequency 
increases monotonically and its derivative with respect to 
edd blows up as e^ 1 * is approached. On the contrary, the 
frequency of the other two modes decrease and their in- 
clinations fall down abruptly in the neighborhood of egj*. 
In the case of A < 1, the monopole frequency changes its 
behavior, which ceases to be monotonie in edd- It grows 
for small edd, but, as the criticai interaction strength is 
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FIG. 9: (Color Online) Low-lying oscillation frequencies for 
triaxial traps as functions of the trap aspect ratio A^ for 
£dd = 1.2 and different values of \ y /\ x . The frequencies are 
normalized by their respective non-interacting values. The 
curves are marked with the colors blue (dark grey), black, 
and red (light grey) corresponding to the highest, moderate, 
and lowest value, respectively. Figs. [9] a-d) show the mixing 
of the two quadrupole-like modes, which go continuously over 
into each other. In Fig. [9] e) a level-crossing in the cylinder- 
symmetric case becomes apparent through an abrupt permu- 
tation of the quadrupole modes (see Fig. [10] for more details). 



approached, it starts decreasing as end grows. The radiai 
quadrupole and three-dimensional quadrupole frequen- 
cies behave as functions of edd for A < 1 qualitatively 
nearly the sanie as for A > 1. For specific values of A, 
though, they might cease to be monotonically decreasing. 
The characteristic feature here is that, for both A < 1 
and A > 1, the three-dimensional quadrupole frequency 
vanishes at eXf, as a signal of global collapse of the gas. 

How the oscillation frequencies depend on A is shown in 
Fig. |8]for edd = 0.8 (dotted) and edd = 1.2 (continuous). 
For A < 1, the quadrupole frequency is larger in com- 
parison to the non-interacting case, while the contrary 
is true for the monopole frequency. As A eventually be- 
comes larger than 1, the monopole (quadrupole) becomes 
larger (smaller) than in the absence of interactions. Con- 
cerning the radiai quadrupole frequency, it turns out to 
be the most insensitive with respect to the dipolar in- 
teraction and is always smaller in the presence of the 
DDI with a minimum around A ss 1.74. The behavior 
of the three-dimensional modes normalized by their non- 
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FIG. 10: (Color Online) Frequencies of the low-lying oscilla- 
tions in units of uj x in the cylinder-symmetric confìguration 
as functions of A for edd = 1.2. The colors blue (dark grey), 
black, and red (light grey) label the frequencies in decreasing 
value. At A ~ 0.94 a level-crossing takes place between the 
two quadrupole modes. In the absence of the interaction, the 
level-crossing happens precisely at A = 1, as can be seen in 
the inset. This difference explains the steep curves in Fig. [5^) 
for 0.94 < A < 1. 



interacting values agrees qualitatively with both dipolar 
BECs [9J and with dipolar Fermi gases in the collisionless 
regime [36|. 



B. Oscillation frequencies in triaxial traps 

In the most general case, Le., in absence of cylinder 
symmetry, the oscillation modes do not to behave as in- 
dicated in the inset of the Fig. [TJ each of the three modes 
becomes a superposition of in- and out-of-phase oscil- 
lations in ali three spatial directions. For this reason, 
the modes are better characterized by their frequencies 
and these are naturally mixed, even if one looks at the 
cylinder-symmetric limit of the triaxial solutions. We 
thus plot the frequencies according to their values and the 
colors blue (dark grey), black, and red (light grey) corre- 
spond to the highest, moderate, and lowest value, respec- 
tively. We exhibit in Fig. [9] the dependence of these fre- 
quencies on X x for different values of X y /X x and edd = 1-2, 
with the frequencies normalized by their respective non- 
interacting values. The situations X y /X x — 2, 1/2, 3, 1/3 
correspond to Figs. la), b), c) and d), respectively. 
These pictures show explicitly that the two quadrupole- 
like modes, denoted by the colors red (light grey) and 
black, are now mixed. If the cylinder-symmetric situ- 
ation is considered, a level-crossing becomes evident at 
A s» 0.94, shown in Fig. [S]e). In contrast to the bosonic 
case [E(j> the DDI affeets the value of A at which the 
level-crossing takes place. For this reason, instead of a 
discontinuous transition, as for dipolar bosons, there is 
a steep continuous line in both the radiai and the three- 
dimensional quadrupole modes for the trap anisotropy 
range 0.94 < A < 1. 
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The level-crossing in the cylinder-symmetric case can 
be seen more clearly in Fig. I1Q[ where the frequencies are 
plotted in units of lj x for edd = 1-2. The vertical line 
marks the level-crossing, which takes place at A ~ 0.94. 
The inset contains a zoomed picture of the frequencies 
for edd = 0. There, the level-crossing happens at A = 1, 
explaining the steep lines which show up in the spectra 
of the quadrupole modes in Fig. [5] e). This shift of the 
level-crossing can be traced back to the Fock exchange 
interaction, which is absent in dipolar Bose-Einstein con- 
densates. 



VII. TOF EXPANSION 

Time-of-fìight expansion experiments are a key diag- 
nostic tool in the field of ultracold quantum gases. In 
BEC's, for example, the effects of the magnetic DDI were 
observed for the first time in 52 Cr by measuring the time 
dependence of the aspect ratios for two different polar- 
ization directions after release from a triaxial trap fili ]. 
In this section, we explore the corresponding problem for 
a strong dipolar normal Fermi gas. 

Dipolar effects are expected to be observed in polar 
molecules, due to their large electric dipole moment. 
Trapping and cooling these molecules requires a strong 
confinement in the polarization axes, to assure robust- 
ness against collapse. Therefore, the suppression of the 
attractive part of the DDI indicates that this system is 
bcttcr dcscribcd by normal hydrodynamics. Initial esti- 
mates of the relaxation time for polar molecules suggest 
that this reasoning remains valid during the whole TOF 
expansion 38]. Thus, we expect the dynamics of the 
dipolar Fermi gas to be described by the equations 



1 d 2 Rj 



K 2 
3A 



£ddQi (R,K X ,K Z 



(73) 



together with the conditions for number conservation 
(|53"|) and momentum deformation (frJU|) . Notice that equa- 
tion (j?3"]) diffcrs from (|rJÌ~j) only due to the absent terni 
—Ri, which is responsible for the trapping potcntial. In 
the following, we discuss the results obtained by solving 
these equation numerically, using the static values of Sec- 
tion [V] for the initial conditions of the parameters Ri (0) 
and Ki(0) as well as R t {0) = and 2^(0) = 0. 

Until now, only axial symmetric traps were involved in 
experimental investigations of dipolar Fermi gases. Nev- 
ertheless, we have learned from studies of dipolar BECs 
how useful triaxial traps can be, for instance in the con- 
text of time-of-fiight experiments. 

Concerning the momentum space, we obtain in the tri- 
axial case similar results as for the cylinder-symmetric 
|, where the aspect ratio K x /K z becomes asymp- 



onc 



totically unity as a result of locai equilibrium in the ab- 
sence of the trap. The anisotropie aspect ratios R X X X /R Z 
and RyXy/R z are plotted as functions of time in Fig. [11] 
for X x = 3 and X y — 5. In the upper and lower graphs 
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FIG. 11: (Color Online) Aspect ratios in real and momentum 
space as functions of time. The dashed curve corresponds to 
the aspect ratio in momentum space K x /K z while the red 
(light grey) and blue (dark grey) curves correspond to the 
real-space aspect ratios R y \ y /R z and R x \ x /Rz, respectively. 
The trap is characterized by \ x = 3 and \ y = 5. The top 
curves correspond to edd = 1 while the bottom ones are for 
edd = 3.5. 



we have set edd = 1 and edd = 3.5, respectively, and we 
find that both aspect ratios become smaller than 1 in the 
course of time. Also for traps with X x < 1 or X y < 1, an 
inversion of the corresponding aspect ratio takes place, 
but in the opposite direction. Such an inversion is typical 
for the hydrodynamic regime and was already observed 
for a two-component, normal Fermi gas with strong con- 
tact interaction [3^, Hl| . 

The value edd = 3.5, chosen above, corresponds to N — 
4 • IO 4 KRb-molecules, with a dipole moment of d « 0.51 
Debye induced by an applied electric field and with trap 
frequencies characterized by uj z — 2it x 280, A^ = 3, and 
X y — 5. Simple arguments like the ones given in Ref. [38| . 
show that the hydrodynamic character of the expansion 
holds, at least, for Tot <C 42. Given that the trapping 
frequencies can be changed at will over a wide range, 
the prospeets for observing hydrodynamic expansion in 
dipolar gases of heteronuclear molecules out of triaxial 
traps are quite promising. 

A further important quantity of the TOF-analysis of 
dipolar Fermi gases is the asymptotic values of the as- 
pect ratios. After the expansion the gas becomes more 
and more dilute and the interaction becomes less and less 
important, even in the case of long-range ones. Never- 
theless, studying the asymptotics in time of the aspect 
ratios may stili be useful because they are approached 
very fast. This is particularly relevant for strong pancake 
traps, where this happens just a few uJ _1 seconds after re- 
lease of the trap, as is shown in Fig. [TT] Although we are 
aware of the inaccuracy of the hydrodynamic approach 
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FIG. 12: (Color Online) Asymptotic behavior in time of the 
aspect ratios as function of the DDI strength edd after release 
from the trap. The dashed curve corresponds to the aspect ra- 
tio in momentum space K x /K z while the red (light grey) and 
blue (dark grey) curves correspond to the real-space aspect 
ratios RyXy/R z and R x ^x/Rz, respectively. 

for small dipole moments, we plot the long-time aspect 
ratios in Fig. Q2] for the whole e^d range. There, we can 
identify the tendency of the DDI to stretch the gas in the 
direction of the applied held in real space, whereas the 
momentum distribution remains always asymptotically 
spherical. 

We would likc to remark that the results presented here 
are in overall disagreement with those obtained by Sogo 
et al. in Ref. (3(| . Translating their findings into our no- 
tation, the ballistic formalism predicts an inversion of the 
aspect ratio only for A > 1. In addition, the aspect ra- 
tio in real space R X X/R Z asymptotically approaches the 
one in momentum space K x /K z for every value of A and 
edd- This disagreement stems from the difference in na- 
ture of both approaches: while hydrodynamics assumes 
locai equilibrium provided by collisions, the ballistic ap- 
proach relies on the assumption of no interaction during 
the expansion. While the latter might be true for weak 
interactions, the former seems to be more adequate for 
strongly interacting gases like the one made out of KRb 
molecules. 



to study equilibrium properties starting from the aspect 
ratios in real and in momentum space and including the 
stability diagram. Apart from that, the hydrodynamic 
low-lying excitations were investigated and a level Cross- 
ing was found in the spectrum, which corresponds to the 
spherically symmetric limit of a cylinder-symmetric trap. 
In the case of a triaxial trap with an external field along 
one of the axes, momentum oscillations were found to 
be two-dimensional and in-phase, just as for cylinder- 
symmetric conhgurations. In addition, we also consid- 
ered the expansion of the gas after release from the trap, 
by solving the equations of motion in the absence of the 
harmonic trap. We found that the characteristic inver- 
sion of the aspect ratio in the course of time after release 
from the trap is also present for dipolar fermions and 
that the fast experimental development of ultracold het- 
eronuclear KRb-molecules makes them quite promising 
candidates for observing these effects. 

We shall like to remark that the theory presented here 
fills an important empty space in the study of normal, 
strong dipolar Fermi gases. Namely, it is applicable in 
the hydrodynamic regime, where collisions provide lo- 
cai equilibrium. Though the prospects for achieving this 
regime with heteronuclear molecules are quite exciting, 
a further gap remains to be filled: Due to the possibility 
of continuously tuning the interaction strength through 
an applied electric field, a theory capable of interpolat- 
ing between the collisionless theory of Ref. [36[ and the 
present hydrodynamic one might be needed in a certain 
range of the parameter space. 

A couple of applications of the present hydrodynamic 
theory for dipolar Fermi gases could be thought of which 
would be useful to provide a deeper understanding of 
these systems. Studying the scissors mode, for example, 
could shed light on the detection of anisotropie superflu- 
idity. The response of the system to a rotational field 
— OL is a further interesting possibility, which allows to 
investigate the moment of inertia of a dipolar gas relative 
to different axes. 



Vili. CONCLUSION 

We have studied both equilibrium and dynamical prop- 
erties of a normal dipolar Fermi gas in a triaxial harmonic 
trap. Using a convenient ansatz for the Wigner phase- 
space function of a normal Fermi gas at very low temper- 
atures, we were able to derive equations of motion which 
govern the momentum and particle distributions as func- 
tions of time, as the trap is shaken or even turned off. The 
dynamical theory developed here allows, as a special case, 
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